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We study creation of bi- and multipartite continuous variable entanglement in structures of coupled quan- 
tum harmonic oscillators. By adjusting the interaction strengths between nearest neighbors we show how to 
maximize the entanglement production between the arms in a Y-shaped structure where an initial single mode 
squeezed state is created in the first oscillator of the input arm. We also consider the action of the same struc- 
ture as an approximate quantum doner. For a specific time in the system dynamics the last oscillators in the 
output arms can be considered as imperfect copies of the initial state. By increasing the number of arms in 
the structure, multipartite entanglement is obtained, as well as 1 — + M cloning. Finally, we are considering 
configurations that implement the symmetric splitting of an initial entangled state. All calculations are carried 
out within the framework of the rotating wave approximation in quantum optics, and our predictions could be 
tested with current available experimental techniques. 

PACS numbers: 



I. INTRODUCTION 

Generally quantum information processing is considered 
theoretically and implemented experimentally as a temporal 
sequence of quantum gates applied to quantum bits that are 
stationary in space. Many classical information processing 
devices function differently however. A microchip usually 
takes as an input currents that are propagating through the 
chip where they suffer classical logic gates. The emerging 
current represent the outcome of the particular classical opera- 
tion that has been carried out. In the present work we consider 
the quantum information equivalent of this way of processing 
information, ie we propagate quantum information through a 
system of interacting quantum systems such that, in the course 
of its propagation through the system, it suffers a non-trivial 
quantum state transformation. In the following we will fo- 
cus attention on harmonic chains for their ready availability 
of exact solutions in non-trivial settings with many excitations 
but see also recent work in spin-systems in this 

direction QUI 

The aim is to design devices to be integrated in a quantum 
information processor in which it will be essential the avail- 
ability of reliable sources of entanglement as well as its trans- 
portation and unitary manipulation inside the processor. At 
the same time the proposed structures can be considered as 
cloning machines of one initial input state. 

This paper is structured as follows. In Section[n]we present 
the basic physical model and its Hamiltonian as well as the 
utilized tools to measure entanglement and its time evolution 
in the system. In section[II]]we describe the Y-shaped struc- 
ture formed by coupled quantum harmonic oscillators, and fix 
the initial conditions for entanglement creation. We present 
and justify the required tuning of the coupling strengths be- 



tween nearest neigbours that creates maximum entanglement 
between the last oscillators in the output arms, and provide 
the resulting unitary transformation between initial and final 
states. We show different results varying the system param- 
eters and compare them with previous ones and with the up- 
per limit for entanglement generation in the system. Then we 
extend the result to multiarm structures in which multipartite 
entanglement is created. Finally we study X-shaped structures 
which work as entanglement beam splitters. In Section llVl we 
study the action of the same structures as a cloning machines 
in which the final states in the last oscillators of the output 
arms can be considered as imperfect copies of the initial in- 
put state. We provide the fidelity expression in terms of the 
squeezing parameter and also extend the results to multiarm 
structures which creates multiple more imperfect copies. 



H. MODEL AND ANALYTICAL TOOLS 

We begin by introducing the basic technical tools that we 
will be employing in this work. 



A. Hamiltonian 

We begin by briefly describing the model system used in 
Ref. 01, a ID chain of interacting quantum harmonic oscilla- 
tors, and the tools used to describe quantum states and entan- 
glement in it. In the rotating wave approximation (RWA) the 
Hamiltonian of the system is given by 



H = 
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(1) 

We may rewrite this Hamiltonian in terms of position and mo- 
mentum operators and pk as is customary in a large part of 
the quantum information literature. Defining a = (q+ip)/y/2 



2 



and at = (q — ip) / V2 we find 



C. Entanglement of Gaussian States 



H = 



i M 



4 flk +fk + \ (Qk+i - Qkf + | (Pk+i -Pkf 
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where c is the coupling constant between next neighbors. In 
this paper we will use its matrix form 



H = -R 
2 



V 
T 



1 M 

R T = ^^QiVijqj +PiTijpj, (3) 



ij=l 



where R T = (qi, . . . , q n ,Pi, ■ ■ ■ ,Pn), and V, T are the po- 
tential and kinetic matrices respectively Holl , In the rotating 
wave approximation we find V = T, 



B. Gaussian states 

Analysing the entanglement properties of infinite dimen- 
sional systems is generally technically involved unless one re- 
stricts attention to so-called Gaussian states (HE!] , which are 
more easily described in phase space introducing the (Wigner- 
)characteristic function. Using the Weyl operator = R 
for £ € R 2 " we define the characteristic function as 



X,(0=tr[pW C ]. 



(4) 



The state and its characteristic function are related to each 
other according to a Fourier- Weyl relation, 



1 
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Gaussian states are exactly those states for which the char- 
acteristic function \p is a Gaussian function in phase space 
0] and are completely specified by their first and second mo- 
ments, d and 7 



X P (0 = Xp(0)e-* 



(6) 



As the first moments can be always made zero utilizing ap- 
propriate local displacements in phase space, they are not rel- 
evant in the context of questions related to squeezing and en- 
tanglement and will be ignored in the following. The second 
moments can be collected in the real symmetric 2n x 2n co- 
variance matrix 7 defined as 



73,* 
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With this convention, the covariance matrix of the n-mode 
vacuum is 7 = l2 n . 

For a Hamiltonian operator like J3j we find that the covari- 
ance matrix of the ground state is lUoll 



7 



Vtv- 1 ® vVt- 1 



(8) 



In RWA we have T = V and the ground state is given by 
7 = In © %i< which is the same as the ground-state of n 
non-interacting harmonic oscillators. 



As the covariance matrix encodes the complete information 
about the entanglement properties of the a Gaussian state, we 
will use it in order to quantify the amount of entanglement 
between two groups of oscillators. In this work we are going 
to use two entanglement measures. When considering pure 
states, the appropriate measure of entanglement is the entropy 
of entanglement which is defined as the von-Neumann en- 
tropy of the reduced density matrix of one subsystem. In the 
Gaussian states setting this reduced density matrix is again de- 
scribed by a covariance matrix 7^4 and the displacements. As 
displacements represent local unitary basis changes they are 
irrelevant when considering the entanglement properties of a 
Gaussian state. If we denote the symplectic eigenvalues of ja 
describing m harmonic oscillators by fa we find 12011 for the 
entropy of entanglement the expression 
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(9) 

Another measure of entanglement that can also be applied for 
mixed quantum states is the logarithmic negativity H2I IT3L 
fl4L fl5l fl6l fl8ll . This measure is easy to compute, possesses 
an interpretation as a cost function ll ill and its monotonicity 
under local operations and classical communication has been 
proven recently 0. Given two parties, A and B, that consist 
of m and n harmonic oscillators respectively, the logarithmic 
negativity is defined as 

m+n 

N(p) = log 2 H^Hj = - log 2 (min(l, | 7i |)), (10) 

3=1 

where p Tb is the partial transposition of p with respect to sys- 
tem B, and 1 1 . 1 1 1 denotes the trace-norm. The second equality 
is obtained employing the symplectic eigenvalues jj of the co- 
variance matrix 7 Ts corresponding to the partially transposed 
state p TB (see Refs. HGU). 



D. The equations of motion 

As we are interested in the entanglement evolution we 
must now find equations of motion for the covariance ma- 
trix. The dynamics of the covariance matrix under a Hamil- 
tonian quadratic in position and momentum operators can be 
obtained straightforwardly from the Heisenberg equation. 



-X(t)=i[H,X}. 



(11) 



For our time-independent Hamiltonian Eq. Q, this leads to 
the covariance matrix at time t as 



r(t) 

where 



exp 



T 
-V 
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These equations of motion will be numerically integrated and, 
in some cases, analytically be solved in this paper. 



HI. CREATING ENTANGLEMENT BY PROPAGATION 

We begin our considerations by devising structures in 
which entanglement is generated automatically during the 
propagation of quantum information. 



A. Bipartite entanglement 

In Ref. |3| a Y-shaped structure as shown in Fig.^was pro- 
posed. One arm, the input arm, consisting of M m oscillators is 
connected to two further arms, the output arms, each consist- 
ing of M out oscillators. Only nearest neighbour interactions 
are considered and it is assumed that the structure is initially 
in the ground state, i.e., at temperature T = 0. At time t = 
we perturb the first harmonic oscillator in the input arm ex- 
citing it either to a thermal state characterized by covariance 
matrix elements J qiqi — j PlPl = z for some z, or a pure 
squeezed state characterized by covariance matrix elements 

7gi9i = l/7pipi = z 




FIG. 1: Sketch of the Y-shaped structure. A single input arm con- 
sists of A/i n oscillators (here 8) is connected to two output arms each 
consisting of M ou t oscillators (here 5). The arrow indicates how to 
group oscillators in the same position n of each branch to apply basis 
change 1 1 81 . 

We observe that for an initial thermal state excitation no en- 
tanglement is ever found between the ends of the two arms of 
the Y-shape. This can be understood from the observation that 
a thermal state is a mixture of coherent states, i.e., displaced 
vacuum states. Displacements do not contribute to the entan- 
glement in a state so that it is sufficient to consider the vacuum 
state. If the system is initialized in the vacuum state however 
it remains stationary and will evidently not lead to any entan- 
glement as long as we are considering the RWA. Therefore 



an initialization in a thermal state cannot yield entanglement 
either. 

For a squeezed state input on the other hand a consider- 
able amount of entanglement is generated. These two obser- 
vations are resembling closely optical beamsplitters which do 
not create entanglement from thermal state input but can gen- 
erate entanglement from squeezed inputs (see Ref. lll9ll for a 
comprehensive treatment of the entangling capacity of linear 
optical devices). 

It is now natural to consider which nearest neighbour cou- 
pling provides the optimal generation of entanglement. These 
optimal couplings can be found numerically by adjusting ran- 
domly individual couplings, accepting changes that improve 
the amount of entanglement that is being generated in the 
structure at a given time for a given input. The results so 
obtained provide the intuition by which one then arrives at 
an analytical argument. In Ref. ^ it was shown that it is 
possible to obtain perfect transmission in a linear chain of M 
oscillators if the couplings are adjusted in a square root law 

V n ,n+1 = V n+ l, n = -yj n(M — 7l) (15) 

V n>n = 1. (16) 

With this potential matrix an interchange between the first and 
the M-th coordinate occurs by waiting for a time t = ir/c (see 
Ref. I2H1 for an analogous argument in spin chains). 

To obtain the equivalent perfect transmission distributed be- 
tween the arms in the Y-shaped structure, we must modify the 
coupling strengths between the oscillator in the junction and 
the first two oscillators in the arms, VM in ,M in +i, in the follow- 
ing way 

V Mia ,M la +i = ^y/M iD (M-M iD )/2. (17) 

In Fig. |2 these couplings are shown for M; n = 8 and 
Mout = 5. To understand the origin of this choice of cou- 
pling strengths we link it to perfect transmission in a linear 
chain. Indeed, with a basis change between both systems it 
can be proved that the choice of coupling for the Y-shape cor- 
responds to perfect transmission in a linear structure consist- 
ing of the incoming arm and one outgoing arm while the other 
outgoing arm is becoming decoupled (see Fig. [3}. Labelling 
the two arms with A, B we take oscillators by pairs in the 
same position n of each arm, and make the following change 
just in arms (see Fig.Q 

^ = -L(Q^ + QS), qB = l^ { QA_QB^ (18) 

with analogous equations for momentum operators p n . Intro- 
ducing this change in the expression for the system Hamilto- 
nian (0 we obtain that A arm gets the perfect transmission 
coupling strengths Eq. dl 51 , while B arm is decoupled. How 
does this allow us to obtain the optimal couplings for the y- 
shape? First, we consider the transformed y-shape depicted on 
the right hand side of fig. [3] We perform the same procedure 
carried out for the linear chain (see Ref. yfj for an analytical 
proof of optimality) to obtain the optimized couplings for this 
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FIG. 2: Square root law profile of the potential matrix non-diagonal 
elements which produces maximum entanglement between both fi- 
nal oscillators in the output arms of the Y-shaped structure (just one 
output arm interactions are shown). The coupling between the last 
oscillator in the input arm, and the first oscillator in the output arms 
must be divided by v2- Mi n = 8, M out = 5 and c = 0.2. 




FIG. 3: Under the transformation <18> the Y-shape on the left hand 
side is transformed into a linear chain consisting of the input arm and 
one output arm while the other output arm is now decoupled. 



situation. Transforming back to the original coordinates, ie 
the situation on the left hand side of Fig. [5] we then obtain 
the coupling strengths suggested in eqs. < !15ll7t . This then 
implies that the initially squeezed state is, in the transformed 
picture, transmitted perfectly to the end of the outgoing arms. 
Inverting the transformation, which is indeed the same trans- 
formation as that implemented by a 50/50 beamsplitter, we 
then obtain a two-mode squeezed state in the original picture. 
The results in fl9il imply that the transformation correspond- 
ing to a 50/50 beamsplitter is the optimal linear optics entan- 
gling transformation for an initial squeezed state. This implies 
that the procedure as a whole, ie the choice of couplings, is op- 
timal - at least in the Gaussian setting. If the initial quantum 
state of the first oscillator is squeezed the final states in the 
last oscillators in the present maximum entanglement at that 
time, as can be observed in Fig. [4] where logarithmic nega- 
tivity is shown. The only effect of increasing the number of 
oscillators (adjusting couplings accordingly) is to narrow the 
peaks without decreasing the maximum entanglement. This is 
contrasted by the behaviour of a system in which all the cou- 
plings between nearest neighbours are equal. Then the magni- 
tude of the entanglement decreases as the number of oscillator 
increases, and irregular time patterns are obtained jj]. 
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FIG. 4: Time evolution of entanglement between the two last oscil- 
lators in the Y-shaped structure when an initial squeezed state is pre- 
pared in the first oscillator of the input arm. Maximum entanglement 
is obtained for t — tt/c. Solid line corresponds to Mi n = M out = 2 
and dashed line to M- m = M ou t = 20. Dotted and dash-dotted 
lines correspond to RWA Hamiltonian with constant couplings for 
Mi n — M ou t = 2 and Mi n = M out = 10, respectively. Coupling 
constant is c = 0.2 and squeezing parameter is z — 10. 



In Fig.[5]we can see how entanglement propagates through 
the output arms in the Y-shaped structure. For each pair 
of oscillators occupying the same position n in each arm 
we measure entanglement at each time. Tuning couplings 
as in Eqs. dl5H7t makes possible maximum entanglement 
production between the last oscillators in each arm for 
t = tt/c. 




t 



FIG. 5: Same as Fig.|4]showing as well propagation of entanglement 
along the output arms between oscillator in the same position of each 
arm. M in = 2, M out = 20, c = 0.2 and z = 10. 
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Perfect transmission means that any dispersion has been re- 
versed in the system at t = tt/c and all the supplied energy 
is accumulated in the two oscillators, one at each end of the 
outgoing arms A and B, with the remaining oscillators in the 
ground state. Therefore this final joint state with 4x4 covari- 
ance matrix 7,4s = 7 a 75 is pure, which is witnessed by 
the following equality [l i] 



-hAB(t = n/c)a} 2 = 11, 
where a is the 4x4 symplectic matrix given by 

iL 2 - 

-IL2 



a = 



(19) 



(20) 



At this perfect transmission instants the evolution of the in- 
put state toward the final one 71 — > jab is an unitary transfor- 
mation. Its expression can be obtained by first considering the 
covariance matrix of the transformed system in which there 
is perfect transmission between initial state and A arm, while 
B arm is decoupled. In the new coordinates Q, P (we take 
out the operator symbols for the sake of symplicity), the joint 
state of the last two oscillators in the arms A, B 
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(Q B P B ) 

(p b2 ) 



(21) 

must therefore equal the direct sum of the squeezed initial 
state (perfect transmission with A arm) and the ground state 
(decoupled B arm ) 



%1B = Tii = 7i © 7i 



z 
1/z 
10 
1 



(22) 



Comparing the two matrices in Eqs. (12 1 i and \22\ we then 
find 



{P AZ ) = 1/z, (Q B ') = {P B ~) = 1 



(Q ) = z, 
(Q A P A ) = (P A Q A ) = (Q B P B ) = (P 



= 0. (23) 



The covariance matrix of the output oscillators in both arms 
of the original coupled Y-shaped structure 



Tab = 7a © 7b 



(24) 

can then be obtained by applying i ll 8t i. Taking into account the 
symmetry (7^ = 73) to notice that (Q A Q B ) = (P A P B ) = 
(Q A P B + Q B P A ) = 0, and with Eqs. we find 



i a2 ) 


(q A p A ) 
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IAB = 1A ® 7B 



z+1 

1/z+l 

z+1 

I/24 



(25) 



This finally leads to the transformation mapping the initial into 
the final state for the Y-shape given by 

7A (t = tt/c) = 7 B (t = tt/c) = i [ 7l (t = 0) + 1] (26) 

For these pure states we can use entropy of entanglement 
Eq.(|9) as a good measure of entanglement fll lisll . To obtain 
the dependence of fa with the squeezing parameter we need 
to determine the positive eigenvalues of the matrix fioll 



B = -iajA, (27) 

where a is the symplectic matrix (1201 . We introduce Eq. (I26i 
(with 71 given by dl4» in Eq. J27> . and take just the positive 
one of the two real conjugated eigenvalues of B 



1 



fa = -Vz+l/z + 2, 



which can be introduced in (|9} with m 
expression S(z). 



(28) 

1 to obtain the 



It is interesting now to compare the obtained entangle- 
ment to the largest possible entanglement the system can yield 
given the total energy Et(z) that have been supplied through 
squeezing of the input oscillator 



E T {z)= l -{(ql) + (p\)) 



l = i(z + l/z)-l, (29) 



where the — 1 stands for a shift in energy such that the en- 
ergy of coherent states be zero, Et(z = 1) = 0. Therefore 
we want to maximize entropy of entanglement 



Mpa log pa) 



i=i 



A l log Ai 



(30) 



with the following constraints for probability and energy 
conservation 



(31) 



i=l 



tr(pH) = ^X i 2E i =E T (z). (32) 

i=l 

The A, are the Schmidt coefficients of the density matrix 
P = Si=i I z) (i I (where we have taken into account sym- 
metry between particles, i.e. = \ib) = and E- L = i 
are the dimensionless eigenstates of one harmonic oscillator. 
Using Lagrange multipliers we obtain the maximum entangle- 
ment the system can yield in terms of the squeezing parameter 
z 



Smax(PA) 



1 

in2 



E T {z) 



In 



E T {z) 



In 



E T {z) 



(33) 



where the total energy is given by Eq. \29\ . Fig.[6]shows the 
obtained entanglement (maximum in time evolution) in the Y- 
shaped structure, Eq. (|9}, and the maximum one, Eq. d33l . in 
terms of z. Both quantities are independent on the coupling 
constant c and on the number of oscillators. 
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FIG. 6: (Solid line) Entanglement obtained in the Y-shaped structure 
between the oscillators in the extremes of the output arms (maximum 
at t — 7r/c) and (dashed line) maximum entanglement in the system 
given a energy E(z), vs. squeezing parameter z. 



B. Multipartite entanglement 

The arguments presented above can be extended easily to 
a structure in which one has Na output arms structure as de- 
picted in Fig.0(see ||5| for a treatment of similar structures). 




FIG. 7: Multiarm structure for N A = 4, M in = 8 and M out = 5. 
Arrows indicates how to group oscillators in the same position n of 
each branch to apply DFT basis change <35> . 



Qr 




~k 2wijk/N A 
Vn ) 



(35) 



with an analogous expression for momentum operators. We 
can now follow the same procedure that yielded Eq. J26l > to 
obtain the unitary transformation which links initial squeezed 
state in the input oscillator with the final state at t = ir/c of 
the last oscillator in one of the Na arms 



lA {t = tt/c) = ^[ 7 i(t = 0) + (N A - 1)1)]. 



(36) 



which reduces to Eq. i2bt for Na = 2. The initial state is 
distributed between Na final states and we compute multipar- 
tite entanglement as the bipartite entanglement between the 
oscillator in the extreme of one arm and the other oscillators 
in the rest of arms. Following the same procedure to obtain 
Eq. d28l we can get now the positive symplectic eigenvalue of 
the reduce covariance matrix in the multiarm structure 

1 



Lh(Na) = -V(z + N A - l)(l/z + N A - 1), (37) 
Na 

which it is introduced in Eq. (|9j to obtain the expected be- 
haviour shown in Fig. |8j multipartite entanglement decreases 
with increasing number of arms. 
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FIG. 8: Multipartite entanglement obtained in the multiarm structure 
(maximum in entanglement evolution at t = tt/c) vs. the number of 
arms Na- 



Now the coupling in the junction, Eq. 1171 . must be 



V Min ,M in+ i = -y/n(M-n)/N A , (34) 

and the basis change which decouples all arms except one 
(for which perfect transmission is obtained) is now a discrete 
Fourier transform (DFT). We apply it to position and mo- 
mentum fi\ operators of oscillators occupying the same posi- 
tion n in each arm (see Fig.0 to obtain the new basis opera- 



C. Beamsplitters 

Extending arguments of previous section we now study X- 
shaped arrays. Such structures may be used to split the en- 
tanglement in a state. In qubit systems, entanglement splitting 
cannot be achieved without loss E3l but we will see that the 
Gaussian setting permits lossless splitting. 

To obtain the equivalent perfect transmission, coupling 
strengths must be adjusted analogously to the Y-shaped case, 
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Eqs. ( I15II16> . but with an interaction applied to the coupling 
between the central oscillator and its four nearest neighbours 
given by 

VM ta -i,Mi. = ^M ln ,M ln +i = ^M- m {M-M- m ). (38) 

As initial condition we now create an entangled two-mode 
squeezed state between the first two oscillators in the input 
arm, labelled with 1 , 2 

7 9 191 = 79292 = 7pipi = 7 P 2P2 = cosh(r) 

79192 = -7piP2 = sinh(r), (39) 

with cosh(r) = z. The resulting time-space evolution of en- 
tanglement is shown in Figs.^|[^ and we sketch an expla- 
nation for the observed patterns in the caption of Fig. |9] 




FIG. 9: a) Under the transformation Eq. 1 1 81 the X-shape on the left 
hand side is transformed into a linear chain and two decoupled arms. 
The initial two mode squeezed state (represented by the yellow disk 
across two arms) is then transformed into a tensor product of two 
single mode squeezed states (individual yellow disks) . One of these 
can move through all the structure, while the other one remains in 
the input arm, giving rise to the complex time-space pattern shown 
in Fig. 1101 b) When there is just one oscillator in the input arm, one 
of the states always remains in the first oscillator and entanglement 
is always present in the input arm as shown in Fig. ll II 



IV. CLONING 

At the times when perfect transmission is obtained adjust- 
ing the couplings, the states at the end of both arms can be 
considered as approximate copies of the input states. How- 
ever, we must be careful to adjust the diagonal element of the 
V and T matrices to optimize quality of the clones as the lo- 
cal rotation of the individual oscillators must now be taken 
into account. The reason is that we are now comparing actual 
quantum states at different times, p±(t = 0) and pA(t = tt/c), 
while in entanglement calculation the local rotations did not 




FIG. 10: Time-space evolution of entanglement in the X-shaped 
structure. For times when one of the squeezed states is in the out- 
put arm, regular transmission patterns are obtained because the other 
one must stay in the input arm. (see Fig. However when the 
state returns to the input arm irregular patterns are obtained as now 
both states are in the input arms. Maximum entanglement is obtained 
when both squeezed states are in the input oscillator. 




FIG. 11: Time-space evolution of entanglement in the X-shaped 
structure when there is just one oscillatr in the input arm (see Fig. 
|9j)). There is always entanglement between the two oscillators in the 
input arm as one of the squeezed states remains always in it. 



affect the amount of entanglement. We find that in Eq. (I16t 
we must modify diagonal elements in the following way 

_ f c/2 for M in + M out even 
Vn ' n ~\ c for M- m + M ont odd ' (W} 
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which can be expressed as 



Vn, 



1 



7[3"(-l) 



(M in +M out )i 



(41) 



It is worthwhile pointing out that with these and any other 
arbitrary value for the diagonal elements we also obtain op- 
timum entanglement creation, as internal degrees of freedom 
are equal for botharms and henceforth entanglement is not af- 
fected. 

The fidelity between the pure input state p\ and one of the 
two copies, described by the reduced density operator pa, is 
given by 



F = = tr(piPA)- 



(42) 



We now insert in this expression the input state p\ = 
\^>i)(^/i\ given by Eq. (|5jl, and consider no displacement 
(d = 0) and initial state Xp(0) equal to one in Eq. 0. Us- 
ing expression @ and solving the resulting Gaussian integral 
we obtain 



^/det(7i + 7I)' 



(43) 



where 71 as the initial covariance matrix of the input oscillator 
(I14t . and 7,4 is the final one at t = tt/c of the A-arm output 
oscillator (7,4 = jb for symmetry). 
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FIG. 12: Time evolution of fidelity between initial squeezed state in 
the first oscillator in the input arm and final state in the last oscillator 
of the output arm in the Y-shaped structure. Maximum fidelity is 
obtained for t = tt/c. Solid line corresponds to Mi n = M ou t = 2 
and dashed line to M in = A/ out = 20. c = 0.2 and z = 10. 



J43t and express 71 elements in terms of the squeezing param- 
eter z, Eq.^J we get 



F(z) 



v/32 + 3/z + IO' 



(44) 



which is depicted in Fig. [O] For z = 1 we have a ground state 
as input and a ground state as output and so the maximum is 
reached at F(l) = 1 and fidelity decreases as squeezing is 
increased. 
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FIG. 13: Fidelity obtained in the Y-shaped structure between the first 
oscillator in the input arm and the last oscillator in the extreme of one 
of the output arms (maximum at t = tt/c), vs. squeezing parameter 
z. 



In the Na multiarm structure the fidelity can be calculated 
inserting the unitary transformation d36l > in the fidelity expres- 
sion i43i to get 



F(z,N A ) 



2N A 



l)(z + l/z) + 2N\ + 2 



(45) 



which reduces to Eq. J44-l > for Na = 2. In Fig. [2|we can 
see that with increasing number of arms we obtain a faster 
decreasing of fidelity than the one obtained for entanglement. 
For Na — * 00 fidelity approaches the value F w 0.575 which 
corresponds to an initial state with 2 = 10 and a final ground 
state. 



V. CONCLUSIONS 



In Fig.^l me fidelity time evolution is shown. The value 
F w 0.575 corresponds to overlap between input squeezed 
state 7 gi9l = l/7 PlPl = z and ground state ja = 1- At 
the same time for which we obtained perfect transmission and 
maximum entanglement, t = tt/c, we obtain as well maxi- 
mum fidelity. At this times rest of oscillators are at ground 
states and relation d26i between covariance matrices of input 
and output oscillators holds. If we insert this relation in Eq. 



We have considered situations in which we propagate quan- 
tum information through a system of interacting quantum sys- 
tems such that, in the course of this propagation, it suffers 
a non-trivial quantum state transformation. In particular we 
have considered entanglement creation, quantum cloning and 
entanglement splitting in systems of quantum harmonic os- 
cillators both analytically and numerically under the rotating 
wave approximation (RWA) as it is appropriate in a quantum 
optical setting. 
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FIG. 14: Fidelity in a 1 to Na doner obtained in the multiarm struc- 
ture (maximum in fidelity evolution at t — n/c) vs. the number of 
arms Na, for z — 10. 



We have presented geometrical structures that allow to cre- 
ate, by propagation, the complex unitary operations required 
for the above processes. Such hard wired pre-fabricated struc- 
tures may also be 'programmable' by external actions to add 
further functionality to them. Therefore quantum information 



would be manipulated through its propagation in these devices 
somewhat analogous to modern micro-chips and as opposed 
the most presently suggested implementations of quantum in- 
formation processing where stationary quantum bits are ma- 
nipulated by a sequence of external interventions such as laser 
pulses. 

All these investigations were deliberately left at a device in- 
dependent level. It should nevertheless be noted that there are 
many possible realizations of the above phenomena. These 
include nano-mechanical oscillators 0], arrays of coupled 
atom-cavity systems, photonic crystals, and many other re- 
alizations of weakly coupled harmonic systems, potentially 
even vibrational modes of molecules in molecular quantum 
computing |23]. We hope that these ideas may lead to the 
development of novel ways for the implementation of quan- 
tum information processing in which the quantum information 
is manipulated by flowing through pre-fabricated circuits that 
can be manipulated from outside. 
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